Zero-bias anomaly in two-dimensional electron layers and multiwall nanotubes 
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The zero-bias anomaly in the dependence of the tunneling density of states u{e) on the energy e of 
the tunneling particle for two- and one-dimensional multilayered structures is studied. We show that 
for a ballistic two-dimensional (2D) system the first order interaction correction to DOS due to the 
plasmon excitations studied by Khveshchenko and Reizer is partly compensated by the contribution 
of electron-hole pairs which is twice as small and has the opposite sign. For multilayered systems the 
total correction to the density of states near the Fermi energy has the form Su/i/o = max(jej, e*)/4eF, 
where e* is the plasmon energy gap of the multilayered 2D system. In the case of one-dimensional 
conductors we study multiwall nanotubes with the elastic mean free path exceeding the radius of 
the nanotube. The dependence of the tunneling density of states energy, temperature and on the 
number of shells is found. 
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I. INTRODUCTION 

Electron-electron interaction results in a singular sup- 
pression of the tunneling (single-particle) density of 
states at. the Fermi surface of low-dimensional metallic 
systemsE The effect, known as the zero-bias -anomaly, 
was first discussed by Altshuler and AronovB for dif- 
fusive systems with a short, range interaction, and by 
Altshuler, Aronov and LccO for the Coulomb interac- 
tion. In the two-dimensional (2D) case the correction, 
to the density of states (DOS) is double-logarithmic,cl 
6v/vq In (|e|r) In (A/e), where r is the impurity scat- 
tering time, e is the energy of a tunneling electron, count- 
ing from the Fermi energy, and A is a high-frequency cut- 
off related to the Coulomb screiening, l/T<cA<g;eF,eF 
being the Fermi energy. ZuzinQ showed that the second 
logarithm in this formula is cut off at low energies for the 
experimental setup of a 2D electron plane screened by a 
metal shield. Rudin, Aleiner and Glazmantl generalized 
the theory of zero-bias anomaly to incorporate the ballis- 
tic energies e > 1 /r, and argued that the correction actu- 
ally has thpj form, 8v jv^ ^ — In^ (ei?/|e|). Khveshchenko 
and ReizerQ analyzed the contribution of collective elec- 
tron excitations, 2D plasmons, to the tunneling DOS and 
obtained an additional correction bvjvQ = (|e| — ei7')/2ei?, 
less singular near the Fermi surface but dominant in the 
wide range of energies. 

In the present paper we consider the interaction cor- 
rection to the tunneling DOS of a different system, con- 
sisting of a number of identical periodically spaced two- 
dimensional electron layers. We also neglect the possi- 
bility of electron tunneling between the layers. While 
analyzing such a setup we primarily have in mind high- 
Tc materials, which are attracting considerable interest 
with respect to the properties of electron-electron inter- 
actions. The electron transport in these materials is ex- 
tremely anisotropic, and the inter-layeg tunneling ampli- 
tude in some crystals can be as weak asQO.05 — 2 i^. Even 



though the motion of electrons tunneling from the tip of 
the tunneling electron microscope (TEM) is confined to 
the outermost layer only, the presence of internal layers 
is important as they participate in the screening of the 
Coulomb interaction. 

We also consider a similar though different tunnel- 
ing geometry which is realized in multi-wall carbon 
nanotubes (MWNT). A typical MWNT consists of a 
few graphite monolayer sheets rolled concentrically into 
cylinders with radius i? ~ 10 nm. At zero doping they 
can be either metals or semiconductors, depending on the 
helical arrangement of the carbon hexagons. In the mea- 
surements of the tunneling DOS the jtyLnncling current 
propagates through the outermost shellElQ while the inter- 
shell tunneling is suppressed. Depending on the degree of 
disorder the transport around the elastic mean free path 
I can be shorter or longer than the radius of the nanotube 
corresponding to the diffusive (I < R) or ballistic {R > I) 
motion along the circumference of the nanotube. For en- 
ergies e exceeding the inverse time of propagation around 
the circumference the zero-bias anomaly is described by 
2D formulas. At lower energies a crossover to the regime 
of a quasi-one dimensional (quasi-lD) conductor is real- 
ized. 

The zero-bias anomaly in diffusive quasi-lD conduc- 
tors was first considered by Altshuler and AronovB. The 
correction to DOS in the lowest order of perturbation 
theory was found to be more singular than in the 2D 
case, Sv/iyQ ^ — l/y^jepr. Working beyond the per- 
turbation approximation, Nazarovll3 found that close to 
the Fermi surface DOS has an exponential behavior, 
ln(z//i'o) ^ |e|~j-(this result was later obtained by Lev- 
itov and ShytovO in a different way). BaJliaiic ID con- 
ductors were studied by several authorsll3E3 and were 
shown to have a power law behavior i'/i'q ^ |e|" of the 
tunneling DOS. The crossover between diffusive and bal- 
listic regimes as well as the temperature behavior of DOS 
in multi-wall carbon nanotubes were recently studied in 
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Ref. |IJ under the assumption that electrons reside on the 
outermost shell only. 

Here we study the zero-bias anomaly due to dynami- 
cally screened inter- and intra-shell Coulomb interaction 
in two- and one-dimensional layered systems and its ef- 
fects on the tunneling DOS. In Section ^ the system of 
two-dimensional layers is analyzed, with both regimes of 
diffusive and ballistic in-plane electron motion consid- 
ered. In Section III we discuss the zero-bias anomaly in 
multi-wall nanotubes assuming that the doping electrons 
are distributed uniformly across the shells. 



II. SYSTEM OF TWO-DIMENSIONAL 
ELECTRON LAYERS 

Let us consider a semi-infinite system of identical con- 
ducting two-dimensional layers separated by the distance 
d (see the Figure in Appendix). The tunneling elec- 
tron from TEM propagates within the upper layer, and 
the inter-layer tunneling is neglected. The properties 
of an isotropic two-dimensional electron system are de- 
scribed by the in-plane Fermi velocity v and the electron- 
impurity scattering rate 1/r. The presence of internal 
layers is important as they contribute to the screening of 
Coulomb interaction. 

The first order perturbation correction to the tunnel- 
ing DOS (sec Refs. of the 2D conductor at zero 
temperature has the form 



I'D 



dLuV{u!), 



kl 



V(w) = 3 



qdq {uj + i/T)U{uj,q)T'^{uj,q) 



(2.1) 
(2.2) 



where is the thermodynamic two-dimensional density 
of states vq = m/ir, counting both spin directions. The 
electron-impurity vertex function is given by 



T-\iU,q) = l- 



i/r 



[{u) + ijrY — <fv 



2„.211/2 ■ 



(2.3) 



The function U{uj,q) stands for the Coulomb interac- 
tion of two electrons residing in the outermost plane and 
screened by the infinite number of conducting layers. To 
find this function, we consider the Coulomb potential 
(j){uj, r) created by the tunneling electron located in the 
outermost plane z = 0. It satisfies the Poisson equation, 
that has the following form in the Fourier representation 
with respect to the in-plane coordinates. 



z — nd) (j){LLj, q, z) 

= AT:e5{z), (2.4) 



n=0 



where the last term in the brackets describes the polar- 
ization charge induced in the system of 2D layers, and 
the polarization operator of a single 2D electron layer is, 

J,, . u^ + ^/r-[{LO + ^/Tf-q'v^Y/^ 

[(^ + ,/,)2_,2,2]l/2_,/, ■ (2-5) 

For the solution of Eq. (^^) we refer the reader to the 
Appendix. We obtain. 



U{uJ,q) = -e(f>{uj,q,0) 



Ane"^ sinh qd 



where k is given by the solution of the equation 
cosh kd ~ cosh qd — K{LU,q) sinh qd, 



(2.6) 



(2.7) 



having a non-negative real part Re fc > 0. Here K(uj, q) — 
2'Ke'^H{u},q)lq is the factor describing the dynamical 
screening of an external charge by a single met allic layer. 
When the wavelength decreases qd ^ oo, Eq. (2.6) gives 



UiLO,q) 



27re' 



9- 27re2n(tj,q)' 



(2.8) 



and the conventional expression for the screened inter- 
action in a single 2D layer is recovered. For static in- 
teractions Eq. ( p^ ) gives U{uj 0,g) = /{q + k), 
where k = is the inverse static screening length. 

In what follows we assume that kc? 3> 1. This condition 
ensures that different layers are (at least for low frequen- 
cies) weakly coupled. In the opposite limit Kd <C 1 the 
system could be treated effectively as the 3d metal with 
the cylindric Fermi surface. This regime is beyond the 
scope of our paper. 



A. Ballistic motion 

First we consider the ballistic limit when the tunnel- 
ing bias exceeds the scattering rate e3> l/r. The main 
contribution to the integral in Eq. {2J1) comes from two 
regions: 1) high-frequency plasmon region, qv lo; 2) 
particle- hole continuum qv > to. 

1) Using the high frequency {qv <ti lo) approxima- 
tion for the polarization operator Eq. ( |2.5| ), H(cj,g) = 
I'aq'^v'^ /2il'^ , where we introduced the notation H.^ — 
Lu(uj + i/r), and solving equation (2.7) with respect to 



Vi(c^) 
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, fdq 2,, 

J q 



1) 



VL^ - uo+u^ - J{Q^ - Lul){n^ -L^-)] ■ (2.9) 



Here uj±{q) represents the upper/lower boundaries of 
the plasmon continuum (of an infinite system of lay- 
ers z = 0,±1,±2,...) that can be parametrized by the 
wavenumber q± in the z— direction. 
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4(9) 



qv 



2 fcoth^i 



1 



tanh ■ 



(2.10) 



The mode i^+{q) corresponds to the uniform charge dis- 
tribution across all layers, q± = 0, while the other mode, 
Lj^{q), describes the the alternating charge in adjacent 
planes {q±^ = n/d). The former mode has a frequency 
gap e* — w(K/(i)^/^, and at gd ^ gives the usual three- 
dimensional plasmon (in anisotropic metal). At — + cxd, 
both branches tend to the usual plasmon spectrum of 
a two-dimensional electron gas. This is illustrated by 
Fig. 1. 




FIG. 1. The plasmon spectrum of an infinite number of 2D 
metallic layers is shown. Dashed line shows the plasmon spec- 
trum of a single 2D layer. Plasmons of different layers interact 
with each other thus creating a band (inner area between solid 
lines). The upper solid line lj = u;+{q) represents the upper 
boundary u) — ui+{q) of the plasmon continuum, with the uni- 
form charge distribution across the layers {q± = 0) while the 
lower solid line marks the lower boundary corresponding to 
the alternating charge in adjacent layers {q± — n/d). 



We chose the momentum cutoff q^^ in Eq. (2.9) to be 
larger than the characteristic momenta of plasmon modes 
but still less than the momentum of particle-hole excita- 
tions, i.e. q^ <C uj/v. 

When the frequency is high > e*, the integral over 
the momentum is dominated by large values of g > \/d 
where uj+ and lo- converge exponentially, lo 
exp(— gd)]. Simple calculations give. 



Vi{uj) = In — 



2) High transferred momenta q ^ uj/v ^ q^^ present 
the contribution of particle-hole pairs. Under the con- 
dition K(i 3> 1, the second term in the right-hand side 



of Eq. (2.7) always exceeds the first term (overscreened 
interaction) and the interaction takes the form. 



U{uj,q) 



1 ^{uj + i/tY 



q^v 



i/t 



y^(a7+Vrp^-<Z^r^ — u — i/t 



(2.12) 



Substituting Eqs. (|2Tg) and (^) into Eq. (gj) and in- 
tegrating from q^ to infinity we find 



V2{UJ) = 



1 



In . (2.13) 



4ei? 2neF(^T q^v 

By ad ding this expression to the plasmon contribution 
(^.11 ) we obtain the following expression for the total 
spectral weight function, 



V(w) = Vi(t^) + V2(CJ) 



1 



1 



1 * 
■In — , w > e*. 



(2.14) 



The results ( 2.11 ) and ( 2.14 ) hold as long as the fre- 
quency exceeds the plasmon gap e*. In this large fre- 
quency and large wave number qd > 1 region the plas- 
mons have virtually no dispersion along the z-direction 
due to weak interaction of electron densities on different 
layers. At lower frequencies and wave numbers, lo < t* 
and qd < 1, electron densities on different layers in- 
teract strongly, and the plasmon spectrum acquires a 
significant dispersion along the z-direction. In this re- 
gion we may approximate = ndq^v^ /A, cu'^ — Kv'^/d, 
exp (2qd) = 1 -f 2qd, and take the inverse interlayer dis- 
tance as the momentum cutoff ^ 2/d. Taking the 



integral in Eq. (2.2) explicitly and extracting the imagi- 
nary part, we find instead of Eq. ( 2.11 ), 



Vi(l^) 



1 

4eF 



(2.15) 



The high-momentum (particle-hole) contribution does 
not depend (as long as Kd ^ 1) on the distan ce be tween 
layers. This can be readily seen from Eq. (2.12) that 
does not contain d. This is quite natural as the particle- 
hole pairs do not induce long-range oscillations of ele ctric 
field. Therefore, we can take the old expression ( 2.13 ) but 
with the new cutoff q^j ^ 2/d, and obtain 



VH = - 



1 



2'KepUJT 



■In — , uj<e*. 



(2.16) 



Integrating the expressions (2.14) and (2.16) with the fre- 
quency according to Eq. (2.1), we found the perturbation 
correction to the tunneling DOS in the ballistic regime 



(2.11) 5vb{e) _ max(|e|,e*) - ep 



1 



4e_F 
In 



AirepT max(|e|,e* 



-] ln[^max(|e|,6*)]. (2.17) 



We observe that the correction (2.17) is less singular near 
the Fermi surface than in the case of a single 2D layer. 

Let us point out that in the limit d — > oo our reault 
differs from the formula of Khveshchenko and Reizer,El as 



the leading contribution (the first term of Eq. (2.17)) is 
twice as small as their result. The difference is due to 
the contribution from particle-hole excitations, V2(w) in 
Eq. (|1|). 
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B. Diffusive motion 

Let us now address the regime of strong impurity scat- 
tering, e <^ l/r, when the motion of the tunnehng elec- 
tron is diffusive and the polarization operator has the 
form, n(w,g) = —VQDq^/{Dq'^ — iuj), here D = v'^t/2 is 
the diffusion coefficient. 

When uj > e = Dn/d, the main contribution to the 



integral (2.2) comes from the region of transferred mo- 
menta UJ / Dk ^ q < y^uj/D, while for ui < e the principal 
contribution is due to yjuj / Dnd < q < y'uj/D. In this 
region the interaction takes the form 



UiLO,q) 



•^oDq^' 



and we obtain 



■In 



Att'^i'qDu! max(cLi, e) 



(2.18) 



For a sufhciently dirty system er < 1. In this case in- 
tegrating over the frequency we obtain for the diffusive 
DOS correction at |e| > e. 



1 



In {D'KS/\e\)ln{\e\T) 



(2.19) 



where the first term is the well-known double logarithmic 
correction of Altshuler, Aronov, and LeeEl and the second 
term is the constant contribution of b allistic frequencies 
w > l/r discussed in the Section [I A. 



For |e| < e, we have instead of Eq. ( 2.19 ) 



1 



ln{Kd) ln(|e|T) + 



In^ er 



(2.20) 



and the first logarithm in Eq. ( ^.19| ) is cut off. For a rela- 
tively clean system with er > 1 (i.e. such that P > d/n) 
the second term in the brackets in this equation is,absent. 

Comparing Eq. (2.20) with the result by ZuzinQ for the 
DOS correction for a 2D electron system screened by a 
bulk metal, we can observe that in the diffusive regime 
at low energies the role of an infinite set of 2D electron 
layers z = d, 2d, 3d... is equivalent to a bulk metal screen 
located at a distance d/2 from the outermost plane. 



be neglected. We treat them perturbatively, in the lowest 
order in the tunneling Hamiltonian, 



d'^p 



2^ J (2^ 



[^l+i(p)+^l_i(q)]V^.(p), (2.21) 



which conserves the in-plane momentum during tunnel- 
ing. Having in mind mainly applications to the system 
with strong in-plane scattering (like high-Tc materials) 
we focus here on the diffusive transport regime e <C l/r. 
The corrections to the density of states to the lowest or- 
der in the inter-layer tunneling amplitude t are shown in 
Fig. 2. Of the six diagrams drawn here the more impor- 
tant ones are the diagrams e) and f) that contain four 
diffusons each. They correspond to the tunneling correc- 
tions to the screened Coulomb interaction rather than to 
corrections to Green's function of the tunneling electron 
itself, which are given by the diagrams a) - d). 





FIG. 2. The diagrammatic representation of the correc- 
tion to the single-particle density of states from the inter-layer 
tunneling processes to the lowest order of perturbation the- 
ory in the tunneling amplitude t. The electron from the mi- 
croscope probe tunnels into white circles. Solid lines stand 
for the Green's functions, wavy lines denote the dynamically 
screened Coulomb interaction, black dots represent the tun- 
neling matrix elements, and the dashed lines are the impurity 
ladders (diffusons). 



C. Interlayer tunneling 

Throughout the analysis performed above we disre- 
garded completely all possible tunneling transitions be- 
tween different layers. We now consider the corrections 
to the density of states originating from interlayer tunnel- 
ing and establish conditions when such processes could 



To compute these diagrams one has to know the inter- 
action Unm{uJ,q) between electrons residing on arbitrary 
layers n and to. As finding such a general expression 
seems to be quite a cumbersome task for a semi-infinite 
system of layers, we utilize the corresponding expressions 
for the infinite system to obtain a qualitative estimate of 
the effect. It has the form (see Appendix), 
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Uoo{uj,q) = 



= UQo{u),q)e' 
27re^ sinh qd 
q sinh kd 



- \n—m\kd 



(2.22) 



The diagrams e) and f) are computed to yield for the 
tunnehng correction, 

V (, A , ,^ » /■ 9^9 t/go(c^,g) tanhfcrf/2 



27r2 



If the frequency is large to > e , the main contribu- 
tion to the integral in Eq. ( 2.23| ) arrives from the in- 
terval 1/d < q < lo/kD, where we can approximate 
C/qo tanh fcc?/2 ~ /v^q^^ to obtain 



Svt{^) _ eH'^VQT |e| 



(2.24) 



For smaller frequencies cj < e, the leading contribution 
comes from the integral < (7 < 1/d and reads, 



(2.25) 



It is worth noting that the inter-layer correction changes 
its sign around the point |e| ~ e and in fact lead to 
increase in the density of states for |e| > e. Compar- 
ing the expr essions ( 2.24 ) and (2.25) with the above for- 
mulas ( ^.19 ) and ( 2.201 ) we observe that inter-layer tun- 
neling correction is small provided that t ^ |e|/(ei?r). 
This much stronger condition than the one that might 
be naively expected (t <C |e|) is due to the enhancement 
of the electron-electron interaction by the inter-layer tun- 
neling. When this condition is violated the summation 
of a wider class of diagrams with all possible transitions 
of diagrams is necessary. 



III. MULTIWALL NANOTUBES 

In this section we consider the interaction suppression 
of tunneling DOS in multiwall carbon nanotubes. A mul- 
tiwall nanotube is built up of M concentric carbon tubes 
(or shells) each of which can be obtained by rolling a 
graphite sheet into a cylinder. 

In tunneling experiments the tunneling current is be- 
lieved to propagate along the outermost shell only while 
the tunneling between the shells is largely suppressed. 
Thus the role of electrons in the inner shells is reduced to 
merely screening the interaction between the outer shell 
electrons. 

The electron band structure of a single metallic car- 
bon nanotube consists of iV (this number is usually of 
the order of a few tens and depends on the doping level 
and gate voltage) conducting sub-bands e„ (A:) character- 
ized by the Fermi velocities along the tube axis w„ and 
around the circumference vx_n- Each of the M shells has, 
in general, its own band structure en(fc)- The electrons 



are scattered between different bands (but mostly within 
the same tube) by impurities, lattice imperfections and 
incommensurate potential of neighboring tubes. Not all 
the tubes are necessarily metallic, some of them could 
well be insulating. Experimental evidence of the internal 
structure of MWNT is usually not easy to obtain and 
leaves some room for speculation. In our previous pa- 
perEj we studied the scenario of the dopants residing in 
the outermost shell only, thus neglecting possible inter- 
action with electrons from internal shells. In the present 
paper we study a different aspect of the problem concen- 
trating on the effects of the finite number M of conduct- 
ing shells on the tunneling DOS of a MWNT. To simphfy 
the problem we consider an approximation in which all 
M shells have the same band structure and the same 
doping level. 

Unlike the case of the 2D conductors considered above, 
in the one-dimensional wires the contribution of plasmon 
frequencies co ^ qv, logarithmically exceeds that of the 
electron- hole interval lu < qv. Therefore the weight func- 
tion can be written in a simple form (compare this with 



Eqs. (2.2) and (2.3)), 



dq 



27r 



^Uoo{uj,q,qm), 



(3.1) 



where the sum is taken over the quantized transverse 
momentum qm — m/Ro, m = 0, ±1,±2... The function 
C/oo(i^: 9, <Zm) should be understood as the 00 element of 
the matrix Uij (w, q) which represents the dynamically 
screened Coulomb interaction of an electron in shell i 
with an electron in shell j and satisfies the matrix equa- 
tion 



C/(c^, q) = F(q) + l/(q) \1{lo, q) [/(c^, q). 



(3.2) 



where ny (a;,g, g^) is the polarization operator which, 
according to our assumptions, is proportional to the unit 
matrix 



n,j {lo, q, q„i) = (5y n = 6ijVi 



(3.3) 



Lu{uj + i/r) 

here we introduced the one-dimensional density of states 
I'l — X]ri(^^n)~^ the average squares of the longitu- 



dinal, and the transverse, electron velocities. 



„,2 _ J2n '"n ^^in 



In Eq. (3.2) Vij{q,m) denotes the bare Coulomb poten- 
tial, 

2e2 r 

Vtj {q, q,n) = / d4)Ko (qRi^) cos mcf). (3.4) 

In this equation we introduced the notations Rfi^{(f>) = 
Rf + Rl — 2RiRk COS0, with Ri being the radii of con- 
centric shells forming the MWNT (i = is the external 
shell, i = M — \ the innermost shell). 
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To proceed further with Eqs. (3.2-^]j) we assume that 
the radius of the j-th shell is a linear function of its num- 
ber Ri = i?o(l — jS,)- In the long-range Hmit qRo <C 1, 



we obtain from Eq. (3.4) for the bare interaction, 



Vtj {q, q„ 



|m| 

In 



(t 



4 max(i,j) 



|m| 



C mm(i,i) 

2^ min(i, j), 



m ^ 0, 
m = 0. 



(3.5) 



Here (3 = 4e~^'' ~ 1.26, with C be ing the Euler constant. 
While the second line in Eq. (3.5) holds provided ^ <C 1, 
the first line utilizes only the approximation qRo <C 1. 

We use the set of eigenvectors w^'^'' (q) an d eigenvalues 
Vk{q, m) of the bare interaction matrix (3^) to write the 
screened Coulomb interaction in the form 



U^jiuj,q,qra) = — T 



wl {m)Wj (m) 



(3.6) 



eigenvalues at characteristic plasmon momenta, q^ 

Lu{LJ + i/T)/Nv^, 



(3.7) 



where we defined the average Fermi velocity v = 
X^n ''^n/N and the dimensionless coupling strength g = 
/(2ttv) in a single channel. The function /Af(e) is given 
by the 00-element of the matrix V^^'^, 



M R^\uj{uj + i/T-^)\ 

the last term representing the contribution of the M — 1 
sound-like plasmons, with values of 7m given in the table. 



The eigenvalues Vfc(g, qm) determine the spectrum of col- 
lective plasmon excitations ujm\q) through the poles of 



the interaction Uij{Lo,q,qm) in Eq. ( |3.6| ), and the eigen- 
vectors w^''^ (m) determine the distribution of charge be- 
tween different shells in these plasmon oscillations. In 
the limit qRo S> 1 one can approximate sum with the 
integral a nd r ecover the usual two-dimensional Coulomb 
potential (2 



Typically, MWNTs exhibit ballistic transport around 
the circumference I = vt ^ 10 — 100 nm > R. At ener- 
gies e below the Thouless energy v/R for the transport 
around the circumference the contribution of the m ^ 
terms is non-singular due to the gap in their spectrum 
and is only weakly dependent on e. The m — plas- 
mons, on th e other hand, are gapless. Their contribution 
to Eq. ( |3.l| ) continues to depend on e and has a singular- 
ity at e — !■ 0. Therefore, concentrating on the low energy 
e < v/R dependence of DOS we can neglect the non- 
singular contribution of the m ^ plasmons and retain 
only the m = term in Eq. ( |3.l[ ). 

The index k labels plasmon modes and equals the 
number of nodes in the charge distribution in the fc— th 
plasmon across the section of a multi-shell nanotube. 
The mode fc = is characterized by the uniform 
distribution of the oscillating charge and corresponds 
to the logarithmic eigenvalue, Vb(g) — e^Mln^^^^, 

with wp"* ~ \/\/M. All other M — \ modes cor- 
respond to g-independent eigenvalues of the bare in- 
teraction, Vk — 2£^e'^ak and therefore have sound-like 
spectrum. The coefficients ak are to be computed 
numerically. They range roughly from a few tenths 
to a few units, e.g. for M = 10 we obtain ak = 
10.0; 2.6; 1.2; 0.72; 0.50; 0.38; 0.31; 0.28; 0.26. The compo- 

(k) 

nents Wq of the eigenvectors are also comp uted nu- 



merically. As the eigenvalues of the matrix (3.6) de- 
pend (at most) only logarithmically o n th e momentum 
g, we can perform the integral in Eq. (3.1) taking these 



M 


1 


2 


3 


5 


10 


IM 





0.35 


0.60 


0.84 


1.38 



The contribution of the fc = plasmon decreases with 
the number of tubes M due to the screening by inter- 
nal shells, while the contribution of sound- like plasmons 
increases roughly linearly with M . Since the numb er o f 
shells may not exceed 1/^ the second term in Eq. ( |3.8| ) 
is never greater than the first one, and at most becomes 
comparable to it at M « 1/^. The approximations > qv 
used above assumes that the case of strong interaction is 
realized, ^ak > 1/^e^vo ~ 1/47V. In the opposite limit 
^afc < 1/4A'', the contribution of the fc— th plasmon mode 
is weak and may be treated per turb atively. 

Substituting the expression ( |3.7[ ) into Eq. (^) and 
evaluating the frequency integral we obtain the follow- 
ing perturbation correction to the tunneling DOS of an 
M- shell mult i- wall nanotube, 



5vM{e) 




/M(k|) + /Af(r-i)ln- 



(3.9) 



We now discuss how the results obtained for multi- 
wall nanotubes could be extended to the very vicini ty o f 
the Fermi surface where the perturbation correction (3.£) 
diverges and the nonperturbative approach is required. 
This question was addressed for a single-wall metallic 
nanotube in Ref. |lj that utili^d the phase approxima- 
tion firat proposed by Nazarovt3 and further developed 
in Ref.Ej. Here we concentrate on the dependence of the 
DOS on the number of shells M skipping the derivation 
that could be found in Ref. |lj. The result for the ballistic 
regime, |e| > 1/r m ainly utilizes the plain exponentiated 
perturbative result (|3.9D, 
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^{e) (X exp 



— In — 

N ^ R\e\ 



/M Re 



(3.10) 



The diffusive regime of lower energies, |e| < l/r, is more 
tricky. We obtain, 



i^(e) oc exp 



(3.11) 



where eo = TTg/2NT. This result holds provided that 
the argument of an exponent is large and the tempera- 
ture small enough T < e^/ep. For higher temperatures 
e^/eo < T < e the tunneling DOS acquires the form. 



i^(r)cxexp<;-1.07j^ 




lni/2 . 



- 7a/ 



ye 



(3.12) 



The crossover between the low-temperature (3.11) and 
the high-temperature ( 3.12] ) regimes can be shown to de- 
pend on the energy through the dimensionless variable 
\e\/y/€oT only. We omit the derivation of the correspond- 
ing scaling function which is quite straightforward and 
follows the same route as in the case when logarithms in 
the exponents are screened off by the presence of a metal 
shield, see Ref. 



IV. DISCUSSION 

We considered the zero-bias anomaly in in the tun- 
neling density of states in layered 2D and quasi-lD ma- 
terials with dynamically screened Coulomb interaction. 
The theory presented above is applicable to high-Tc ma- 
terials and semiconductor heterostructures as well as to 
the multi-wall carbon nanotubes. We showed that the 
presence of many conducting shells in 2D systems weak- 
ens the singularity in the DOS correction at the Fermi 
level: for energies below e* ~ nD/d in diffusive systems 
or below the plasmon gap ^(K/d)^/^ in ballistic systems 
the DOS correction becomes logarithmic oc ln|e|, rather 
than double-logarithmic cx In^ |e| as in the case of a single 
layer. 

In addition to the logarithmic terms mentioned above 
that originate from the broad region of the transferred 
frequency and momentum there are also terms com- 
ing from the regions around the plasmon and particle- 
hole singularities. These contribute correction to DOS 
max(|e|, e*)/4ei? that is not singular but dominant in the 
wide range of energies |e| > ln(ei?r)/T. The contribution 
of the particle-hole continuum is especially interesting in 
2D ballistic systems. It turns out that the contribution 
of the particle-hole continuum to the DOS correction. 



which corresponds to the region uj < qv in the interac- 
tion frequency/momentum transfer, of the same order as 
the contribution originating from the plasmon pole but 
has a different sign. Namely, the particle-hole continuum 
increases the density of states. The total effect is still a 
suppression of DOS for low energies as the plasmon con- 
tribution is always larger than that of the particle-hole 
continuum (twice as large for a single 2D layer). 

In derivation of DOS in 2D systems we assumed that 
the ground state is the free electron model. It is pos- 
sible, however, to incorporate also the Fermi liquid in- 
teraction in the calculations of the tunneling DOS. In 
particular, the Fermi liquid theory predicts that thermo- 
dynamic density of states is renormalized vq Zvq by 
a Z-factor (0 < Z < 1) describing a step in the quasi- 
particle partition function. We obtain however that the 
relative correction 6v/vq = \e\/4,ep due to Coulomb in- 
teraction is in fact universal as it is independent of the 
Fermi liquid interaction strength Z . 

The validity condition for neglecting the intcr-shell 
tunneling was found for the diffusive regime. It was 
shown that instead of the naively expected requirement 
of the tunneling time being greater than the charac- 
teristic measurement time , the much more stringent 
condition ^ epT/e 3> e^^ must be satisfied. 

With respect to multi-wall carbon nanotubes our anal- 
ysis was concentrated on the dependence of DOS on the 
number M of coaxial shells for both the regimes of per- 
turbative and strong suppression. We found that this 
dependence can be described by a single function 7m 
found numerically. The energy and temperature depen- 
dence of DOS found in Ref. |l^for a nanotube with a sin- 
gle conducting shell is preserved for an arbitrary number 
of shells. 
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APPENDIX A: SOLUTION OF THE POISSON 
EQUATION 



To solve the Poisson equation (2.5) we first derive the 
Green function for the infinite problem, 
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47re^n 



OO 
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n— — OO 



d{z-nd)j G{z,z')^d{z-z'), 

(Al) 



with the boundary condition G{z, z') —f when 



00. The homogeneous solutions of Eq. (Al) are easily 
written in terms of quasiperiodic Bloch functions, 



V'± (z) = e**^"'^ (sinh [q{z -{n + l)d]- e^^'^ 

X sinh - nd)]) , nd<z<{n^\)d (A2) 

with fc given by the solution of the equation 

27re^n 

cosh kd — cosh qd sinh qd (A3) 

having a positive real part, Re fc > 0. This choice means 
that the function 'ip-{z) decreases and the function ip+{z) 
increases with z increasing. The Green function can then 
be written in the conventional form 



G{z,z') 



1 



i}-{z)il;+{z'), z>z', 

2q sinh qd sinh kd \ ip+{z)ip^{z'), z < z' . 

(A4) 



To find the solution for the semi-infinite problem, 
Eq. ( ^.4| ), we impose a fictitious charge Q located at 
z = -d/2 (see Fig. 3). 



z=-2d 



z=-d 
z=-d/2 

+ z=0 

z=d 
z=2d 
z=3d 



FIG. 3. The crystal consisting of a semi-infinite system of 
2D metallic layers {z = 0,d,2d,..) is shown. The thin solid 
line ai z — —d/2 hosts a fictitious (image) charge Q that 
accounts for the absent z = — d, —2d, .. layers. 

This charge has to be found from the condition that 
the total potential of the electronic charge — e and the 
fictitious charge Q 



contains only the exponent e"^^ for the negative values of 
the coordinate, —d/2 < z < 0. 

This condition ensures that the electric field in the 
outside region {z < 0) decays exponentially with the 
distance from the outermost plane, with the fictitious 
charge, therefore, taking care of the absent {z = —dn, 
n = 1,2, ...) planes. We obtain the following value of the 
fictitious charge 



V'-(o) 



r^kd 



^+{-d/2) e 



-kd , 



(A6) 



Substituting E q. (|A6 | ) into Eq. (A_5) and making use o f 
Eqs. (A2) and ([A4|) we obtain the final expression ( |2.6| ). 



In Section II C the formula for interaction Unm of two 



electrons residing at different layers z = nd and z = md 
in an infinite system is used. To obtain it one can write 
Unm = ~4:Tr e'^G{nd,md), which leads to the expression 
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